Clustering transitions in vibro-fluidized magnetized granular materials 
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We study the effects of long range interactions on the phases observed in cohesive granular materials. 
At high vibration amplitudes, a gas of magnetized particles is observed with velocity distributions 
similar to non-magnetized particles. Below a transition temperature compact clusters are observed 
to form and coexist with single particles. The cluster growth rate is consistent with a classical 
nucleation process. However, the temperature of the particles in the clusters is significantly lower 
than the surrounding gas, indicating a breakdown of equipartition. If the system is quenched to 
low temperatures, a meta-stable network of connected chains self-assemble due to the anisotropic 
nature of magnetic interactions between particles. 



Vibro-fluidized inelastic spherical particles are an im- 
portant model system which capture the essence of dis- 
sipative interactions on the statistical properties of dry 
non-cohesive granular materials [|l|. This system has 
emerged as an important test bed to investigate the ap- 
plicability of dissipative kinetic theory [^,|[ . Experiments 
measuring the position and velocity of individual parti- 
cles show the formation of clusters, non-Gaussian veloc- 
ity distributions (to varying extent), and the violation 
of equipartition [||j5|. In a number of applications, addi- 
tional cohesive interactions often exist due to the pres- 
ence of moisture, electrostatic screening, and magnetiza- 
tion §. 

In this Letter, we introduce a novel system consisting of 
magnetized particles inside a vibrated container. This en- 
ables us to study the effect of long-range attractive inter- 
actions on the formation of clusters and the velocity dis- 
tributions of cohesive granular materials. Non- magnetic 
particles are also present in our system to define a system 
temperature that depends on the vibration amplitude of 
the container. A gas-like phase of magnetized particles 
is observed at high system temperature. If the system is 
slowly cooled below a transition temperature, compact 
clusters precipitate and grow in time. If the tempera- 
ture is rapidly quenched, an extended mesh of particles 
is observed to self-assemble due to the anisotropy of the 
interaction. The structure and growth of the clusters is 
consistent with both a classical nucleation process, and 
recent theories of the dipolar hard sphere model [0-0. 
However, the velocity distributions and the associated 
granular temperature of the free and clustered particles 
show the influence of dissipation and cohesion. 

The idealized interaction between two dipolar hard 
spheres separated by distance r is defined as 
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U = Uhs + ^{fii ' iTj) - ^{fTi ■ fij){lTj • fij), (1) 

where Uhs corresponds to the hard core repulsion inter- 
action, /i is the dipole moment, and r is the interparticle 
vector connecting the centers of dipoles i, j. This type of 
interaction has been put forth as an extremum model for 
ferro-fluids [^,12|. If the potential is averaged over all 



possible dipole arrangements, the familiar van der 
Waals interaction is recovered. This implies that such 
systems should have a well defined liquid gas transition. 
However, extensive simulation and theoretical work on 
the dipolar hard sphere model have not conclusively de- 
termined the existence of critical phase transition. Re- 
cent experiments on ferro-colloid suspensions suggest the 
existence of a liquid gas transition, but open questions 
still remain Other experiments using nickel coated 
glass micro-spheres in a liquid solution have found the 
existence of ring and mesh-like states [|4j. 

The particles used in our experiments are chrome steel 
spheres with a diameter of a = 0.3 cm (with a high degree 
of sphericity Sa/a ^ 10~^) and mass m = 0.12 g. Each 
sphere has been placed in a ramped field of 1 x 10^ G to 
embed a permanent moment of /i ~ 10~^ emu per par- 
ticle. The apparatus consists of a 30.0 cm diameter flat, 
anodized aluminum plate, with 1.0 cm sidewalls and a 
clear acrylic lid, that is weakly coupled through a rigid 
linear bearing to an electromechanical shaker. The sys- 
tem is leveled to within 0.001 cm to ensure that the plate 
is uniformly accelerated. The measured acceleration of 
the plate F = Auj'^/g, where A^uj are the amplitude and 
angular frequency and g is the acceleration due to grav- 
ity, is varied between F = — 3.0 ^, at a; = 377 rad s~^. 
Image data is acquired through a high speed Kodak SR- 
1000 digital camera with a spatial (temporal) resolution 
of 512 X 480 pixels (250 f.p.s.). The covering fraction 
of magnetic particles, 0, defined as the ratio of the area 
of the particles to that of a close packed mono-layer, is 
varied from cj) = 0.01 0.15. We also place glass 
particles of equivalent mass, with a fixed volume fraction 
of (j)p = 0.15 to introduce additional stochasticity into 
the system. The glass particles act as thermal carriers 
to ensure that as the magnetic particles condense, the 
temperature of the system does not go to zero. 

The granular temperature is defined as Tg = ^m(v^). 
By measuring velocity components along the two hor- 
izontal axes Vx^y we build a distribution of veloci- 
ties whose width, given by the second moment, deter- 
mines the granular temperature. The form of the v^^^^- 
distributions deviate from a Gaussian consistent with 
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FIG. 1. Examples of chains (a), rings (b), and crystallites 
(c) observed. Rings appear be the most stable configuration. 
The scale bar denotes 1 cm. (d) Snapshot of the system at 
(j) — 0.09, where T is lowered to Ts from the gas state after 
1092 s. (c) The system at = 0.15 after a rapid quench from 
the gas state into the network state. Images (c,d) are the size 
of the system. 



previous observations |^]. The system temperature T, is 
determined by the velocity distribution of the gas of both 
magnetic and non-magnetic particles. We will return to 
a more detailed analysis of the velocity distributions for 
the gas and clustered phases after discussing the nature 
of the clusters. 

Our experiments begin by maintaining a gas of mag- 
netic and glass particles at a particular (j). By lowering T, 
a fraction of the magnetic particles begin to self-assemble 
into the simple structures shown in Fig. |l|(a-c). We de- 
fine the temperature at which this occurs as the tran- 
sition temperature, Tg. These structures are transient, 
and over time will evolve into large clusters that coexist 
in steady state with single particles. An example where 
a single ring has grown in time to form an extended clus- 
ter, is shown in Fig. |l|(d) [|l5[. However, clusters of a very 
different conformation are produced if T is quenched far 
below Tg. Figure |l|(e) is an example of the system after 
a rapid quench. This extended mesh will eventually re- 
arrange to take on more a compact form, therefore we 
consider this phase to be a meta-stable version of the 
clustered phase. 

To characterize the gas and clustered phases observed, 
and the nature of the transitions, we plot the phase dia- 
gram [Fig. ^ |l^. The connected line shows the transi- 
tion temperature Tg as a function of (j). Above Tg, short 
lived dimers and trimers are observed along with single 
particles, but at and below this temperature, small seeds 
of magnetic particles [Fig. ^(a-c)], precipitate from the 
gas. Initially seeds take on the form of short chains that 
quickly become either rings or more compact states (it is 
noteworthy to mention that rings > 4 particles are en- 
ergetically more favorable than a chain of equal particle 
number |l|,|l3). Fi gure |l|(d) shows an image of the sys- 
tem with 6 = 0.09 after 1092 s at Tg. If T is increased 
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FIG. 2. (a) The phase diagram of temperature T versus 
the covering fraction of the particles (f). The plate acceler- 
ation r, scaled by gravity, is also shown for clarity. A gas 
phase consisting of single particles and short lived dimers and 
trimers are observed above a transition temperature Tg that 
depends on 0, shown by the solid points. To evaporate a 
cluster in the gas phase one must go past Ts denoted by the 
hysteresis region. Below Ts, dimers and trimers act as seeds 
to the formation of compact clusters that coexist with single 
particles [Fig. 0(d)]. If T is rapidly quenched from the gas re- 
gion to very low T highly ramified networks of particles form 
[Fig. Se)]. 



above T^, clusters show hysteretic behavior. This hys- 
teresis depends on the ramping rate of T, and is not ob- 
served to disappear over laboratory time-scales. There- 
fore it appears that the observed phase transition is first 
order. 

In the gas phase, we measure both the velocities of the 
non-magnetized glass particles and the magnetized steel 
particles and find that their distributions are indistin- 
guishable. Figure || shows the distributions of velocities 
at ^ = 0.09 for each particle species, viz. the gas of glass 
particles, the gas of magnetized steel particles (prior to 
and after nucleation of a cluster), and the particles within 
a cluster. The data shown corresponds to the system in 
Fig. 0(d), although the results are the same for differ- 
ent values of (j). When particles nucleate into clusters 
the fluctuations in the particle positions compared to the 
surrounding gas dramatically decreases, as shown by the 
narrow distribution in Fig. |^. However, the surrounding 
particles still have a temperature identical to that of the 
pure gas phase. The ratio of the system temperature to 
the cluster temperature Tc is T/Tc ~ 60. This effect is 
inherent to the non-equilibrium nature of this system, 
where effects due to friction limit the thermal agitations 
expected in the constituent particles that comprise gels 
and polymers [0. The existence of a finite measurable 
temperature in the clustered phase does imply however 
that particles are still mobilized, which can be seen by 
the substantial re-arrangements that occur during cluster 
growth. 

By utilizing the Hoshen-Kopelmen flq] cluster identi- 
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FIG. 3. The probabihty distribution functions P(v) of the 
velocity components in the gas and clustered phases on a log 
linear scale, (o) The p.d.f.'s for the magnetic particles in both 
the gas and clustered phases. The narrow distribution is from 
particles within a cluster and the broad distribution (of the 
same symbol) is for the surrounding gas. (□) Magnetized, 
and (•) glass particles in the gas phase. The velocities are 
absolute velocities given in (cm s~^) not rescaled by the tem- 
perature. The solid line is a Gaussian fit to the magnetized 
gas distribution. 



fication algorithm individual clusters that form at Tg are 
identified. We are then able to track each clusters' size, 
center of mass, and radius of gyration. This method was 
also employed to measure the temperature of a single 
cluster and the surrounding gas. The average number 
of particles contained in a cluster (ric) = nc/n, where n 
is the total number of clusters formed, is plotted versus 
the radius of gyration Rg [Fig. ^. n is unity for most 
times, but may vary at very short times. This data was 
obtained for 1092 s immediately after T was lowered to 
Tg. Using the scaling relation ric ~ Rg we mea- 

sure the compactness of the clusters. Figure ^ clearly 
shows a crossover in the dimensionality of the clusters. 
Clusters with Rg < 5.0 <j follow closely to the spatial di- 
mension. As the number of particles per cluster increases 
with time the trend for the fractal dimension, at all ^, 
is more consistent with (i = 1.4 ± 0.1. This implies that 
clusters are initially more compact and become more ex- 
tended as they grow. The change in dimensionality can 
be understood by observing the growth of a cluster. As 
particles join the cluster they do so by forming highly 
mobile chains that are tethered by one end to the sur- 
face of the cluster. If the free end of the chain is able 
to connect back onto the cluster an excluded area is es- 
tablished. During this process free particles may become 
trapped by the chains, which in time can produce rear- 
rangements and re-opening of these regions. This can be 
seen by a close inspection of Fig. |l|(d) and in Ref. ]T5[ . 

Next we plot the cluster growth rates by measuring the 
number of free particles at Tg for various (j). We find that 
the average radius Tc given by nl^^ (n = 1 for Fig. || over 
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FIG. 4. Average cluster size {ric) versus the radius of gy- 
ration Rg/cF. The scaling (ric) = aR^^ where a = tt^ 
shows a crossover ai Rg ~ 5 a indicating that the clusters 
start more compact and become more extended. 



all time) contained in a cluster grows linearly as the sur- 
rounding free particles are depleted to a saturation value 
determined by 0. The initial linear growth of Vc is con- 
sistent with the Wilson-Frenkel growth process of hard 
sphere colloidal crystals, where individual particles are 
absorbed onto the surface of a crystallite |20|, although 
the mechanism is quite different. The change in the lin- 
ear growth of Tc to a more arrested rate is caused by two 
effects. As the clusters grow the number of free particles 
is reduced and therefore the growth rate must decrease 
as the system approaches steady state. Data sets shown 
in Fig. ||(a) correspond to clusters that nucleate close to 
the center of the cell. Clusters are mobile and eventually 
drift to the edge of the cell [see Fig 0(d)]. Due to deple- 
tion forces |2l|, the cluster will never be able to return 
to the center of the cell once in contact with the side- 
wall. As a result, half of the surface is not absorbing free 
particles, therefore the growth rate is further decreased. 

To characterize the temperature of each cluster and 
emphasize its difference from the system temperature, 
we measure the cluster temperature Tc as a function of 
time at (/) = 0.09 [Fig. ^(b)]. Tc is not only observed to 
be much lower than the surrounding gas phase, but also 
observed to decrease with increasing cluster size. 

Finally, we discuss the meta-stable network phase 
that is observed within the clustered phase, shown by 
the shaded region of the phase diagram [Fig ^. This 
phase [see Figure 0(e)], directly demonstrates the inher- 
ent anisotropy of the potential described in Eq. This 
phase is attained by quenching to a T far below Tg from 
the gas state. When the energy of the system is dra- 
matically decreased in this manner, the particles quickly 
form chain configurations because locally the head-to- 
tail alignment is energetically most favorable. However, 
as discussed above, when the length of the chain is > 
4 particles, the energy of the chain can be further low- 
ered by bending into a more compact structure. At low 
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FIG. 5. (a) The radial growth rate Tc/ct at T = Ts for 
each (j). The sohd hnes are fits to the data over a time when 
the cluster is not in contact with the wall. The apparent 
crossover in growth rate is due to a finite size effect. Linear 
growth is consistent with a Wilson- Frenkel growth process, 
where individual hard spheres are absorbed onto the cluster. 
Symbols for each = 0.03 (A), 0.05 (□), 0.07 (o), and 0.09 
(•). (b) The cluster temperature Tc as a function of time for 
(j) = 0.09. The temperature of the surrounding gas remains 
at a constant value of T = Ts —2.2 g (cm s""*^)^ 

temperatures the probability of overcoming the potential 
barrier of bending is decreased, therefore we are able to 
observe the meta-stable state that is comprised mostly of 
string like objects for a longer period of time. It is self- 
consistent that the deeper the quench, the longer the net- 
work phase survives, which is compounded even further 
by the dissipative interactions between particles. How- 
ever, it is interesting that this phase shows similarities to 
recent Monte-Carlo simulations of quasi-2D dipolar hard 
spheres [|2|, where dissipation is not considered. 

In summary, we visualize the nucleation of clusters in 
a non-equilibrium system by using magnetized granular 
materials. We have directly investigated the nature of 
the clusters formed by measuring their temperature, di- 
mensionality and their growth rates. The temperature of 
the clusters is found to be considerably lower than the gas 
phase indicating a breakdown of equipartition. We also 
find that the phenomenology of the clusters is similar to 
that shown by recent work on dipolar hard sphere mod- 
els, including a networked phase [|l0|,|2|. Although our 
system is weakly dissipative, it may provide still useful 
insights on the nature of phase-transitions in ideal dipo- 
lar system due to our ability to directly visualize and 
tune the particles in the experiments. 
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